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Abstract. Recently two-loop electroweak corrections to the neutral current four-fermion processes at high 
energies have been presented. The basic ingredient of this calculation is the evaluation of the two-loop 
corrections to the Abelian vector form factor in a spontaneously broken SU(2) gauge model. Whereas 
the final result and the derivation of the four-fermion cross sections from evolution equations have been 
published earlier, the calculation of the form factor from the two-loop Feynman diagrams is presented for 
the first time in this paper. We describe in detail the individual contributions to the form factor and their 
calculation with the help of the expansion by regions method and Mellin-Barnes representations. 



1 Introduction 

Electroweak higher order corrections in the high en- 
ergy Sudakov regime ^|z] have recently attracted a new 
wave of interest 3,4,5,6 J ,8, iCDF 

the upcoming colliders, the LHC and an International Lin- 
ear Collider, for the first time the characteristic energies 
i/s of the partonic processes will be far larger than the 
masses of the W- and Z-bosons, M\y,z- In view of the 
expected experimental accuracy, one has to take into ac- 
count radiative corrections at the two-loop level which are 
enhanced by up to four powers of the large electroweak 
logarithm ln(s/M^ z ). These are present in virtual cor- 
rections to exclusive reactions like electron-positron or 
quark-antiquark annihilation into a pair of fermions or 
gauge bosons. 

For the high energy behaviour of the neutral current 
four-fermion processes the analysis of the leading loga- 
rithms (LL) in [5] was extended to the next-to-leading 
(NLL) and next-to-next-to-leading logarithmic (NNLL) 
level in |31llUlllllll2| . With the help of evolution equations 
which describe the dependence of the amplitude on the en- 
ergy, the logarithmic corrections were resummed to all or- 
ders in perturbation theory in NNLL accuracy. Neglecting 
the fermion masses and the mass difference between the 
W- and Z-boson, the logarithmically enhanced part of the 
two-loop corrections to the total cross section and to vari- 
ous asymmetries was obtained including the ln n (s/M^ v z ) 
terms with n = 4,3,2. The results up to NLL accuracy 
have been confirmed by explicit one-loop [Sill 5111 711 18| and 
two- loop 20,21,23,26 calculations. 

a Bernd Feucht in publications before 2005 



For energies in the TeV region the subleading logarith- 
mic contributions are comparable in size to the leading 
terms due to their large numerical coefficients. Thus the 
calculation of the remaining two-loop linear logarithms is 
necessary to control the convergence of the logarithmic ex- 
pansion. These corrections represent the next-to-next-to- 
next-to-leading logarithmic (N 3 LL) contributions. In con- 
trast to the higher powers of the electroweak logarithm, 
they are sensitive to the details of the gauge boson mass 
generation, in particular they depend on the Higgs boson 
mass M#. Thanks to the evolution equations, the NNLL 
calculation involves only massless Feynman diagrams at 
the two-loop level. But the linear two-loop logarithm re- 
quires the evaluation of vertex corrections with the true 
masses of the gauge bosons and the Higgs boson. 

In [22, 27, 33, 52] the previous analysis is extended to 
N 3 LL accuracy. The application of the evolution equation 
approach to the linear two-loop logarithm and the nec- 
essary ingredients are described in detail in |34j . From 
the viewpoint of loop calculations, the most complicated 
contributions are the massive two-loop corrections to the 
Abelian vector form factor which will be defined in sec- 
tion[21 In |33U34j the form factor results are used together 
with the evolution equations to obtain the N 3 LL two- 
loop corrections to the four-fermion scattering amplitude 
in a spontaneously broken SU(2) gauge model. The ad- 
ditional infrared-divergent electromagnetic contributions 
are separated according to the prescription developed in 
|27j. Finally the effect of the mass difference between the 
two heavy electroweak gauge bosons W and Z is taken 
into account by an expansion around the equal mase case, 
whereas a value of the Higgs mass identical to Mw is suf- 
ficient for the desired accuracy. In this way electroweak 
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corrections to the total cross section, forward-backward 
asymmetries and left-right asymmetries of the neutral 
current four-fermion processes are obtained including all 
large two-loop logarithms and leaving an estimated theo- 
retical uncertainty of a few per mil to one percent for the 
production of light fermions [3"5II54] . 

This calculation and the discussion of the results are 
not repeated here. The following sections are instead ded- 
icated to details of the loop calculations needed for the 
form factor mentioned above. The paper is organized as 
follows. In section[21the Abelian vector form factor in the 
spontaneously broken SU (2) gauge model is denned. Then 
section 01 describes the evaluation of the two-loop ver- 
tex corrections to the form factor. Contributions from the 
renormalization of the fields, the coupling constant and 
the gauge boson mass are added in section 0] Finally we 
discuss the result for the form factor in section [5] and con- 
clude with a summary in sectional The appendices list the 
Feynman rules for the SU (2) model (appendix^) and ex- 
plain two important methods used in our calculation: the 
expansion by regions (appendix[5J) and the Mellin-Barnes 
representation (appendix Q. At last appendix IT)1 lists the 
contributions in a theory with a mass gap which are neces- 
sary for the separation of the electromagnetic corrections. 

2 The Abelian vector form factor 

The Abelian vector form factor F determines the fermion 
scattering in an external Abelian field. It is the factor 
which multiplies the Born term J-£ = l 4 , (Pi)l >J " l P(P2) in the 
corrections to the Abelian vector current = F(Q 2 )J r ^, 
where pi denotes the outgoing and pi the incoming fer- 
mion momentum and Q 2 = — (pi — p-i) 2 - At high energies, 
we consider the Sudakov limit fJH] Q 2 — ► oo, so Q 2 ^> 
M 2 for every gauge boson or Higgs mass M, we neglect 
fermion masses, p\ — p\ = 0, and we omit terms which 
are power-suppressed by at least one factor M 2 /Q 2 . 

The four-fermion amplitude A describes the neutral 
current scattering // — > f'f of a fermion-antifcrmion 
pair into a different fermion-antifermion pair. At high en- 
ergies and fixed angles, where all kinematical invariants 
are of the same order and far larger than the gauge boson 
mass, s ~ |t| ~ \u\ ^> M 2 , it can be decomposed into the 
form factor squared and a reduced amplitude A, 

a = imL f ^a, (i) 

s 

where g is the weak SU(2) coupling. The collinear diver- 
gences appearing in the limit of vanishing gauge boson 
mass, M — > 0, and the corresponding collinear logarithms 
are known to factorize. They are responsible, in partic- 
ular, for the double-logarithmic contribution and depend 
only on the properties of the external on-shell particles, 

I43| . Thus, for each fermion-antifcrmion pair of the four- 
fermion amplitude the collinear logarithms are the same as 
for the form factor F discussed above. The reduced ampli- 
tude A in Q therefore contains only soft logarithms (cor- 
responding to soft divergences in the limit M — > 0) and 



renormalization group logarithms. It can be determined 
with the help of an evolution equation (43U44U45) . 

The decomposition of the four-fermion amplitude A 
and the calculation of the reduced amplitude A are de- 
scribed in detail in piUUllU12j . The NNLL approxima- 
tion of the two-loop contribution to the form factor F 
can be obtained with the help of another evolution equa- 
tion |11U12| . whereas the N 3 LL result including all large 
logarithms requires the two-loop calculations with massive 
gauge bosons presented in the following sections. 

The form factor is calculated as a real function of the 
variable Q 2 > 0, i.e. in the Euclidean region. For its appli- 
cation to the four-fermion amplitude described above, the 
analytic continuation to the Minkowskian region s > 
according to Q 2 = — (s + i0), where iO denotes an in- 
finitesimal positive imaginary part, leads to the substitu- 
tion ln(Q 2 /M 2 ) = \n(s/M 2 ) - m. 

The calculation is performed in a spontaneously bro- 
ken SU{2) gauge model. Reference jSJ discusses in detail 
all effects resulting from the difference of this model with 
respect to the standard model of particle physics, where 
the isospin SU{2) group for left-handed fermions is mixed 
with the hypercharge U{\) group through the mass eigen- 
states of the Z-boson and the photon. 

In contrast to the standard model particles and 
Z, we work with the neutral SU{2) gauge bosons W a , 
a = 1,2,3, which all have the same mass M = Mw- 
The generators of an SU(N) gauge group in the funda- 
mental representation are labelled t a . Their Lie algebra 
involves the structure constants f abc . The Casimir opera- 
tors of the fundamental and the adjoint representation are 
C F = (N 2 - 1)/(2N) and C A = N, respectively. In addi- 
tion Tp — 1/2 is needed. In the special case of an SU(2) 
group the generators t a correspond to half the Pauli ma- 
trices, and f abc = e abc , C F = 3/4, C A = 2. We prefer to 
use the general symbols t a , f abc , Cf, Ca and Tp instead 
of their specific SU(2) values in our calculations. This has 
the advantage that we can easily convert the results to the 
case of the hypercharge i7(l) gauge group. 

The Feynman rules of the vertices needed for our cal- 
culation are listed in appendix [X] We use the Feynman- 
't Hooft gauge, where the masses of Goldstone bosons and 
ghost fields are equal to the gauge boson mass M and the 
gauge boson propagators have the form —ig^/i^k 2 — M 2 ). 
We work with the Lagrangian as a function of the unrenor- 
malized quantities, so instead of calculating diagrams with 
counter terms, we have to replace the bare mass and cou- 
pling constant in the one-loop result by the corresponding 
renormalized quantities as described in section 01 

3 Vertex corrections 

In this section the two-loop vertex corrections to the 
Abelian vector form factor are presented. The Feynman 
diagrams contributing to the vertex corrections are de- 
picted in figures ri][ll and Solid lines with arrows denote 
fermions, wavy lines denote gauge bosons, short-dashed 
lines with arrows denote ghost fields and long-dashed lines 
stand for the Higgs boson (H) or Goldstone bosons (4>), 
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Fig. 1. Fermionic vertex correction 




(a) (b) 




(c) (d) 
Fig. 2. Abelian vertex corrections 

depending on the labels. The graphs in figures Et), Eli) 
and|2K) also appear in a horizontally mirrored version, so 
their contributions have to be counted twice. 

The three figures group the Feynman diagrams in sub- 
sets which are separately gauge-invariant when adding 
the corresponding renormalization contributions from sec- 
tion The fermionic contribution of the graph in figure ^ 
is proportional to n/, the number of fermions running in 
the closed fermion loop. Figure El represents the Abelian 
graphs (in addition to figure Q which are present also in 
an unbroken U(l) theory like QED. Finally figure |3 shows 
the non-Abelian graphs, which include the contributions 
from the Higgs mechanism. The Abelian contribution only 
counts the part of the graphsEt") andEt) which is propor- 
tional to Cp. The other part of these two graphs, which 
is proportional to CfCa, belongs to the non-Abelian con- 
tribution. 

The fermionic contribution has been calculated exactly 
in 22, i.e. for all Q 2 , not only Q 2 > M 2 , showing the 
good agreement of the Sudakov limit with the exact contri- 
bution for energies larger than 300 GeV. The high-energy 
asymptotic limit of this result is quoted in section 13.21 




(e) (f) 
Fig. 3. Non-Abelian vertex corrections 

Note that we state in this paper the individual vertex cor- 
rection, self-energy correction and renormalization terms, 
whereas in |22j only the total fermionic contribution to 
the form factor is given. 

The Abelian graphs have been evaluated in N 4 LL 
approximation, i.e. including all large logarithms and 
the non- logarithmic constant. These calculations are pre- 
sented in sections 13.31 to 13.61 The non-Abelian graphs, 
especially figure [3Ji) , are more complicated to evaluate, 
as they have three massive propagators each (compared 
to two for the Abelian graphs). We have only evaluated 
them in N 3 LL accuracy as the non-logarithmic constant 
is not needed for the insertion of the form factor result 
into the four-fcrmion amplitude. The corresponding cal- 
culations can be found in sections 13 . 71 and 13.81 

3.1 Reduction to scalar integrals 

From each Feynman vertex diagram in the figures ECU 
by applying the Feynman rules in appendix ^ we get a 
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vertex amplitude of the following form: 



■pi 



(2) 



where ip(Pi,2) are doublets of Dirac spinors in the SU(2) 
isospin space corresponding to the incoming and outgoing 
fcrmion, and is a quadratic matrix both in the spinor 
space and in the isospin space (for each Lorentz index /i). 

For vanishing fermion masses (in the Sudakov limit) 
the vertex amplitudes can be written as T% = FvFq, 
where is the Born amplitude and F v is a contribu- 
tion to the form factor. The scalar quantity F v can be 
extracted from the vertex amplitude by projection (see 
e.g. 051): 

p = Trjj^r^) 

AN{d-2) Pl - P2 ' [) 

where dimensional regularization 02] is used with d = 
4 — 2s as the number of space-time dimensions, and N = 2 
for SU{2). The trace runs over the spinor and the isospin 
indices. By applying this projection, we get a linear com- 
bination of scalar loop integrals. For convenience, we sep- 
arate the integration measure as follows: 



li 



-v 



d d k 
7 2^Y 



47T 



M 2 



: (M 2 



d d k 
I^ 2 



Here fj, is the mass scale of dimensional regularization, 
a — g 2 /(47r) with the weak coupling g, je is Euler's con- 
stant and S e = (47r) e e _£7B . Within the MS renormaliza- 
tion scheme, S £ is absorbed into fi 2e by a redefinition of 
/i, and as we set [i = M in the end, the prefactor in front 
of the square brackets gets especially simple. 

The reduction of the Feynman amplitudes to scalar 
integrals has been performed with the computer algebra 
program FORM 0B]> and the evaluation of the scalar in- 
tegrals, as described in the following sections, has been 
done with Mathematica |4*9"] . 

We have not performed a reduction of the scalar inte- 
grals to so-called master integrals by a method like inte- 
gration by parts 50, 5 J, as the number of scalar integrals 
obtained from the Feynman diagrams is not too big and 
most of the scalar integrals can easily be evaluated in a 
semi-automatical way starting from our expressions for 
general powers of the propagators, which are presented in 
the following sections. 



3.2 Fermionic vertex correction 

The fermionic vertex correction of figure ^ h & s been eval- 
uated in |22| . where the integration of the inner fermion 
loop has been done first, leaving a one- loop integral feasi- 
ble by the standard Feynman parametrization technique. 
The contribution of the fermionic vertex correction to the 
form factor is 



/ a 

F v>nf = C F T F n f l^— 



M 2 



2 s 



s 2 



V 20 8 , 29 
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Fig. 4. Scalar graph for planar vertex correction 
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(M 2 



(5) 



where £ = ln(Q 2 /M 2 ), and Cs ~ 1.202057 is a value of 
Riemann's zeta function. 



3.3 Planar vertex correction 

The reduction (see section 13. If) of the planar Feynman 
graph in figure leads to scalar integrals correspond- 
ing to the graph in figure The numbers enumerate the 
inner propagators and correspond to the indices i of the 
propagator powers rii and of the inner momenta fcj, the di- 
rections of which are indicated by the arrows. Solid lines 
stand for massive propagators, dashed lines for massless 
ones. 

Apart from propagators in the denominator, one scalar 
product remains in the numerator which cannot be ex- 
pressed linearly in terms of the denominator. We have 
chosen this irreducible scalar product to be 2k§ ■ (fcg — k§). 
The set of scalar integrals is then covered by the following 
function [k = feg, £ = k§ — kg): 

F hA (n u ...,n 7 ) = e 2e ^ (M 2 ) 2£ (g 2 )"-"- 4 
d d k f d d £ (2fc-£)™ 7 



ijT d/2 J ilT d/2 (p _ 2pi . £)m (£2 _ 2p2 . t)m 
1 

{k 2 - 2pi • k) n * {k 2 - 2p 2 ■ A;)™* (k 2 - M 2 )«= 
1 



((k - I) 2 - M 2 )«e ' 



(G) 



with n = ^123456, where we use riy... = rij + rij + . . . as a 
shorthand notation. The scalar integrals have been defined 
in such a way that they do not carry a mass dimension. 

The evaluation of the scalar integrals has been per- 
formed with the expansion by regions (see appendix IBl . 
The following regions contribute to the planar vertex cor- 
rection 1521: 



(h-h): 
(lc-h): 



Pi, 



Q 
Q 
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(2c-h) 
(lc-lc) 
(2c-2c) 
(h-s') 



k\\p 2 , l~Q 

fell Pi, i\v\ 

fc||P2, l || P2 

k ~ Q, fcg = k — I 



M 



By k ~ Q we mean that each component of the vector fc is 
of the order of Q. And k \\ pi indicates a region where the 
momentum k is collinear to the external momentum pi : 



k || Pi 



M 2 

Q ' 



Q: 



M, (7) 



M 2 



k\\p 2 ^> k + ~Q, k_~—, fc x ~ M, (8) 

where k± = (2pi^ ■ k)/Q denotes the components of k in 
the direction of p 2 and p\ respectively, and the vector k±_ = 
k — (k-/Q)pi — (k + /Q)p2 is made up of the components 
of k perpendicular to pip. 

The leading term in the expansion of the (h-h) re- 
gion corresponds to the massless integral with M = 0, 
which is well known 53,54,55 . The (h-s') region is of or- 
der (M 2 / Q 2 ) 2 ~ n6+e and therefore suppressed by at least 
one factor M 2 /Q 2 with respect to the (h-h) region for all 
scalar integrals we need, i.e. for m < 1, i = 1, . . . , 6. So we 
do not need to consider the (h-s') region. The (lc-lc) re- 
gion is of order (M 2 /Q 2 ) 4 ~" 1356+ ™ 7 > the (2c-2c) region of 
order (A/ 2 /Q 2 ) 4 "" 245l3+n7 . Both are suppressed if n 7 > 0, 
i.e. if the numerator is present, and are only evaluated 
for 717 = 0. The leading contributions from the (lc-h) 
and (lc-lc) regions can be expressed by one- and two-fold 
Mellin-Barnes representations (see appendix IH)l : 



F 



(lc-h) 



LA 



(m,...,m) = [ 



T(f - n 3 )r(f - 7^i6 + n 7 )r(n 35 - f ) 

r(n 1 )r(n2)r(n 3 )r(n 5 )r(n 6 )r(d - n 126 + m) 

ioo 

x J ^-.r(-z)r(i-n 26 -z) 

— ioc 

r(n 6 + z)r(n 37 - n 4 + 2;)r(ni 2 6 — 4 + z) 
v - 

r(| - ra 4 + n 7 + z) 



(9) 



F^ c - lc) (n 1; ...,n 6 ,n 7 = 0) 



1 



Q 2 " 



IOO zoo 



dzi /" dz2 



r(n 1 )r(n 3 )r(n 5 )r(n 6 )r(^-n 2i ) J 2m J 2m 

— zoo —ioo 

r(- Zl )r(n 13 - j- gi)r(| - "i + *i) 

r(f-n 4 + ^i) 

x r(f - n 24 + zi) — - d 

r(| -n 5 -z 2 ) 

x r(| - 7745 - Z2)-T(ni3 56 — d+ z 2 )r(n 5 + z\ + z 2 ) . 



(10) 



For symmetry reasons we get f^A ^ from Fj^'^ and 



^LA 
724. 



(2c-2c 



from F, 



(lc-lc) 



LA 



by exchanging ri\ <-» 772 and 713 



The integration contour of the Mellin-Barnes integrals 
runs from —ioo to +700 in such a way that poles from 
gamma functions of the form _T(. . .+z) lie on the left hand 
side of the contour ("left poles") and poles from gamma 
functions of the form r(. . . — z) lie on the right hand side 
of the contour ("right poles"). 

The Mellin-Barnes integrals in © and (|l()(l are solved 
by closing the integration contours either at positive 
or negative real infinity and summing over the residues 
within the contour. The integrals develop singularities 
at points in the parameter space of the n% where a left 
pole and a right pole glue together in one point. Some of 
these singularities are cancelled by zeros originating from 
gamma functions in the denominator, e.g. in F^ when 
n 6 — 0. Here the result is given by the limit n 6 — > 
to which only the residue of the integrand at z — or 
z = —tiq contributes. 

Other singularities in the parameter space are can- 
celled between several regions. This is the case for the 
pole 1 / (n 3 — 714) which is cancelled between the (lc-h) and 
the (2c-h) regions. Another pole 1/(7113 — 7124) is cancelled 
between the (lc-lc) and the (2c-2c) regions. Such singu- 
larities, which are regularized analytically with the pa- 
rameters rii in individual regions, are typical for collinear 
regions in the Sudakov limit. The sum of the contribu- 
tions from all regions is well-defined in the framework of 
dimensional regularization. 

In some cases, the first Barnes lemma (see appendixlUll 
is used to solve one of the two Mellin-Barnes integrations 
in HlOfl . In more complicated cases first all residues which 
produce singularities are extracted, and the limits of the 
analytic regularization and of dimensional regularization 
(e — > 0) are performed before summing up the remaining 
residues. These sums are then solved by Mathematica 
or looked up in a summation table (e.g. in |56j). 

By adding together the contributions from all regions 
we have obtained the results for all scalar integrals origi- 
nating from the reduction of the planar Feynman diagram. 
As examples, we show the results for the scalar graph with 
all propagators present and various powers of the numer- 
ator: 

F LA (1, 1, 1, 1, 1, 1, 0) = ^C 4 + y C 2 - 6C 3 £ + ^ 4 , 



F L A(l,l,l,l,l,l,lHy£-10C 3 , 

F LA (1, 1, 1, 1, 1, 1, 2) = ^ + i (-£ + l 



12 2 



(ii) 

(12) 
(13) 



Here and for all other results of individual scalar integrals, 
we omit the specification u + 0(e) + 0(M 2 /Q 2 )" of the ne- 
glected terms. The result ljTT|) has already been calculated 
in 
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Fig. 5. Scalar graph for non-planar vertex correction 

The complete Feynman diagram in figureEK) involving 
contributions from all scalar integrals with different rij 
yields the following planar vertex correction: 



(2c'-2c): k 3 || pa, £ || p 2 

(us'-us'): h ~ M 2 jQ, fc 4 ~ M 2 /Q 

(lc-us'): k\\p u k 4 ~M 2 /Q 

(us'-2c): k 3 ~M 2 /Q, £ \\ p 2 

The leading term of the (h-h) region is known from 
the massless case |53U54II55] . As for the planar ver- 
tex correction in the previous section, the (lc-lc) and 
the (2c-2c) regions are of order (M 2 /Q 2 ) 4 ~™ 1356+ " 7 and 
(M 2 /Q 2 ) 4_ ™ 2456+ri7 respectively. They are suppressed for 
n 7 > and are therefore evaluated only for 777 = 0. The 
leading contributions from the regions, apart from (h-h), 
can be written as one-fold Mellin-Barnes integrals or sim- 
pler expressions: 



Fv,la = C% ( 



Utt 



2 f l_ 

M 2 



2f 



-C 2 +3£- |tt 5 



2e 2 
11 

~ T 

-32(3 - 7T 2 



T ' c 



8 4 ™> 13 2 41 

— 7T 4 + 62C 3 H 7T 2 

15 S3 12 8 



'M" 



(14) 



F NP '{^,...,717) = [ — 



M 2 



2 \ 2-Tt 35 +£ 



r(| - 7724)^(1 - 77 16 + 77 7 )J>35 ~ f ) /" dz_ 

r(n 1 )r(n 2 )r(n 3 )r(n 5 )r(d-n 1 246+n 7 ) 2 J 2ni 

— ioo 

r(-z)r(| - 77146 - - "4246 + ^37 - z) 

r^-n w -z) 
r(m + z)r(f - 773 + z)r(77 1246 - f + z) 
r(m 6 + z) 



(16) 



F^ F '(m, ... ,n 6 ,n 7 = 0) = ( 



3.4 Non-planar vertex correction 

The non-planar Feynman graph in figure^) involves the 
scalar integrals depicted in figure[S| With the choice 2k^-ke 
for the irreducible scalar product, the scalar integrals are 
written as 



^ N p(77 1 ,...,77 7 )=e 2 ^ (M 2 f E (Q 2 ) 



— p 2e~t E ( M 2\2e ( n 2^n-m-A 

d d k f d d £ (2k -tf 



l7T d/2 J i7r d/2 _ k _ £)2)m (( p2 _ k _ g)2y l2 
1 



(k 2 - 2pi ■ fc)™3 (p - 2p 2 ■ £) ni (k 2 - M 2 )"= 
1 



(P - Af 2 )™« 



(15) 



with k — fcg and I = k§. The following regions contribute 
to the non-planar vertex correction [581152) : 



r(n 16 - f )r(77 1356 - d) 
r{ ni )r(n 3 )r(n 5 )r(n 6 )r^-n 24 ) J 2m 



d r(jl b - ?7 4 + Z)r(f - 77j + z) 
X I 7713 — % — Z) =-, 

2 r(n 156 -n 4 -l + z) 

T(f - 7724 + Z)ra - 77 34 + 2) 



r(f - 774 + z) 



(17) 



2 \ 4-Tt 3 45 6 



r(lc-2c)^ m _\ — I M \ „-niTT 2ej E 



F^ p >(m,...,n 7 ) = — 



r(n 37 - n 2 )r(n i7 - rei)r(| - ni 3 )r(| - 7724) 

r(77 3 ) J r(77 4 ) J r(77 5 ) J r(776) r(| - n 12 + n 7 ) 2 



x r(77 35 - |)r(77 46 - f ) 



(18) 



-fNP (77l,...,77 7 ) = [ — j 



2 \ 4 — «2345 



g -ni7T e 2e-f E 



(h-h) 


k < 


- Q, 




j Q 


X 


(lc-h) 


k 


pi, 


£r. 


j Q 


(h-2c) 


k > 


- Q, 


*\\ 


P2 




(lc-lc) 


k 


pi, 


t\\ 


Pi 


X 


(2c-2c) 


k 


1 P2, 


^11 


Pi 




(lc-2c) 


k 


Pi, 


e\\ 


P2 


X 


(lc-lc') 


k 


Pi, 




1 Pi 



r(d - 77234)^(7724 - f )r(?72345 - d) 

r(n 2 )r(n 3 )r(n 4 )r(n 5 )r(± - 77 16 + n 7 )r{d — 772347 J 



^(1 - 77 2 + Z)r(| - 77 16 + ?7 7 + Z) 
T(| - 776 +^7 + 2) 



(19) 
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(us'-us') / \ 

'(ni,...,n 7 ) 



M_ 2 

Q 2 



2\ 8— »ii256— 2n 34 — 2e 



X e 



2£7E rjd - n 13i )r(d - n 234 )r(n 13 - |K>24 - |) 

r(ni)r(n 2 )r(n 3 )r(n4)r(n 5 )r(n 6 ) 

x /> 23 45 - d)r(n 1346 - d) , (20) 



F, 



(lc-us') 



("1 



NP V 1 , " • ' ' U \ 

r(| - » 4 )r(f - 7i 13 )r(d - n 234 )/> 24 - f) 
x J r(n 2 )r(n3)r(n4)r(n 5 )r(n 6 )r(n 4 7 - m)r(f - n 3 ) 

x r(n 46 - f )r(n 347 - f )/> 23 45 - d) . (21) 

Using the symmetry of the non-planar graph under the 
exchange of the parameters ni «-> n 2 , n 3 «-> n 4 and 



n 5 «-> n 6 , one gets .F N p 

p(lc-lc) p(2c'-2c) , p 



(h-2c) 



from F 1 , 



(lc-h) 



XP 



4 2 P C " 2C) from 



(lc-lc) 7-,t2c'-2c) c t-^(Ic-Ic') i 7-,fus'-2c) c 7-^flc-us') 

*. i w\ i en i anc [ p\ ) f rom pv 



NP ' NP „, NF NP uum 1 NP 

The expression <|20[) for the (us'-us') region is valid for 
general nr although it does not involve wj explicitly: The 
only dependence on n-j of this region is cancelled by the 



prefactor (Q 



03- 



We have checked the completeness of our set of re- 
gions by writing the full scalar integral for arbitrary pa- 
rameters rii (except n-j = 0) as a four-fold Mellin-Barnes 
representation. From this expression, we have extracted 
the residues yielding the non-suppressed contributions and 
have found 11 terms with exactly the same dependence on 
M 2 /Q 2 as the 11 regions listed above. 

The evaluation of the Mellin-Barnes integrals is done 
as described in the previous section. The structure of sin- 
gularities needing analytic regularization is more compli- 
cated than in the planar case. Various poles involving 
combinations of the parameters rii are cancelled between 
the collinear regions (lc-lc), (2c-2c), (lc-2c), (lc-lc') and 
(2c'-2c). 

The contributions of all regions sum up to the results 
for the scalar integrals originating from the reduction of 
the non-planar Feynman diagram. Examples of these re- 
sults are 

F NP (1, 1, 1, 1, 1, 1, 0) = ^£ 4 - ^C 2 + 20C 3 £ - ^tt 4 , 

(22) 



r 2 _4 
-4 ,r r 2 



F NP (1, 1, 1, 1, 1, 1, 1) = -£ 4 - — C z + 14C 3 £ - — , (23) 
4 6 90 

F NP (1, 1, 1, 1, 1, 1,2) = 4 + -(-4C + 7) + ]c 4 - C 3 

7T 2 \ 7T 4 7T 2 

— +9 £ 2 + (14C 3 - 30)£ - — -4C 3 ' 



!)() 



38. 



(24) 



^(1,1,1,1,1,1,3) = ^ + V-7£ " ' 



W - ^ 3 

4 2 



7T 2 59 

IT + T 



C 2 + [ 14C 3 - — I £ 



211 



it 4 _ 3 2 571 

9o- 6C3+ r + - 



(25) 




Fig. 6. Scalar Mercedes-Benz graph 



The result (12211 for the scalar graph without numerator 
is known from |58| . The complete non-planar vertex cor- 
rection with contributions from all scalar integrals is as 
follows: 



F 



v,NP 



= c 



2 \ 2e 



-n 2 + 12)C 2 



-£ 4 - -C 3 
3 3 



40( 3 + - 28 J r - ^tt 4 - 72C 3 



2si +0{e) + o(^r 



(26) 



3.5 Vertex correction with Mercedes-Benz graph 

Figure illustrates the scalar integrals resulting from the 
reduction of the Mercedes-Benz graph in figure |2t). With 
our choice of 2p 2 ■ k§ as the irreducible scalar product, the 
scalar integrals are defined as 



F BE (ni, . . . ,n 7 ) = e 
d d k f d d £ 



2s lE ( M 2)2e (g2y>-, l7 -4 

(2p 2 • k)" 7 



in d/2 J ilT d/2 (£2 _ 2pi . l)m {p _ 2p 2 . £)"2 
1 

X (k 2 - 2pi • k - M 2 )"3 (p - M 2 ) n * (k 2 ) n * ((jfc - e) 2 ) n » ' 

(27) 

with k = &5 and I = k 4 . The list of relevant regions for 
the Mercedes-Benz graph is shown here: 



(h-h) 


k < 


-Q, 


I' 




(lc-h) 


k 


Pi, 


t' 


-Q 


(h-2c) 


k < 


- Q, 


I 


P2 


(us-2c) 


k < 


- M 2 /Q, 


I 


P2 


(lc-lc) 


k 


Pi, 


I 


Pi 


(2c-2c) 


k 


Pi, 


i 


P2 


(lc-2c) 


k 


Pi, 


I 


P2 



The leading contributions of all regions could be evaluated 
for general n-j, but the (us-2c) and the (2c-2c) region are 
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only non-suppressed for n 7 — 0: 



ii,i 2 ,i 3 >0 

r{ni + i 3 )r(n 4 + i 2 ) 



l'l23<«7 



i\\i 2 \ h\ (n 7 - U23)! r(ni)r(n 2 )r(n 3 )r(n 4: ) 

r(% - n 35 )r(d - ni3456 + ii)r(n 12 34, 56 ~ d) 
r(n 5 )r(n 6 )r(d - n 356 + ii2 3 )r(^d - ^123456 + n 7 ) 



dz 
2tt7 



• / —r{-z)r{i -n h& + i 12 - z)r{n h + z) 



x r{& - n 2 4 + n 7 - ii2 + z) 



r(n 3567 - i 123 - I + z) 

("13567 - «12 - 2 + Z ) 

r(d - n 36 - 2n 5 + zi 23 - z) 
- n 36 - 2n 5 + zi2 - z) 



^BE" h3 ( n i'---> n 7) = ( : ^ 2 - 



Af 2 



-ni7r 2£7e 



x r(n 35 - f )r(f - n 146 ) 
r(ni)r{n 2 )r{n 3 )r(n & )r(d - n 124& ) 

ioo 

x I ^-r(-z)r(l-n 2A6 -z) 

— ioo 

r{n 6 + z)r(f - n 5 + n 7 + z)r(n 1246 - j + z) 
r(|+n 7 + z) 



(29) 



be (,ni,...,n 7 J - [ -^2 



e -ni7r g2E7£ 



r(| - n 2 )r(| - n 35 )r(| - n 1356 + n 2 ) 
r(n 2 )r(n 3 )r(n 4: )r(n 6 ) 



-T(f - n 56 + n 7 )r(n 24 - f )/>356 - f ) 
- n 356 + n 7 ) 
f2 \ 6— T12346 — 2715+117 — £ 



r(d - ni 356 )/ n ((i - n 3 56 + ra 7 ) 



(30) 



77i(uS-2c) / . / 

be '{ni,...,n 7 ) = [ — 



AP 



r(| - n 5 + n 7 )r(d - n 256 + n 7 )r(n 35 - f ) 



r(n 2 )r(n 3 )r(n4)r(n 5 )r(n e )r(n 5 - Tin) 

x r(n 25 - n 17 - I )r(n 56 - n 7 - f )r(n 245 6 - n 7 - d) , 

(31) 

/ Jl,/-2 \ 4—7113456 



^be (ni,...,n 7 ) 



ii2 <nj 



n 7 \ 



r(ni + i 2 ) 



i\\i 2 \ (n 7 - ii 2 )\ r(m)r(n 3 )r(n 4 )r(n5) 



zoo too 

dzi f dz 2 



r(n 6 )r{!± -n 2 + n 7 ) J 2m J 2ni 



n-zi) 



—ioo — too 



r(n m -h- zi)-T(ni356 - h - d - zi) 

r(m + «2 - zi) 



x r(n 4 + i x + zi)r{-z 2 )r{!± - n 56 + n - z 2 ) 
r(n 1 - n 5 +i 2 - z 1 - z 2 ) 



x r{n 5 + z 2 ) ■ 



r(ni - n 5 - z x - z 2 ) 



x r(| - ni + n 7 - « 2 + 21 + 22) 
r(| - n 2 + rt 7 + zi + z 2 ) 



T(| + n 7 + Zl + z 2 ) ' 



(32) 



; (ni,...,n 7 ) = ( — 



2 \ 4-n 2 456+"7 



Zl23<«7 

X V 



r{ni + ii)r(n 37 - ij) 

. «i!«2!i 3 ! ("7 - «123)! r(ni)r(n 2 )r(n 3 ) 

^(n 2 - n 137 + i 2 ) r{± - n 35 + h)r(§ - n 6 + ^23) 

r(n 4 )r(n 5 )r(n 6 ) r(l - n 13 )r(d - n 356 + i 123 ) 



x r(d - n 256 + n 7 + i 3 )r(n 56 - i 23 - f ) 



(28) x r(n 24 56 -n 7 -d) 



j^flc^c) , \ / "' 

^be [n 1 ,...,n 7 ) = [ — 



(33) 



M 2 
2 



g-niw g2£7-E 



r(n 2 -ni 6 )r(7|-n2)r(7 



^56 



n 7 )r(n 24 - |) 



r(n 2 )/> 3 )/> 4 )r(f - n 16 )r(| - n 6 + n 7 ) 
x r(n 35 - I ) . (34) 

For the summation indices we use the shorthand notation 
*12— — *i+«2 + - ■ an d the multiple summation is defined 
in the following way: 



112- .. <"7 



2l,22,--->0 Zl— 22=0 

We were able to reproduce the above expressions H28JI - 
(I34|l for the regions by writing the full scalar integral with 
general TLi clS Si triple sum over a three-fold Mellin-Barnes 
integral and extracting all non-suppressed contributions. 
Our evaluation of the (h-h) region is in agreement with 
the known results for the massless diagram |53II54II55| . 

The (lc-lc) region is of order (M 2 /Q 2 ) -1 when n\ — 
fta = 714 = 715 = Jig = 1. But in all these cases the inverse 
power of M 2 is cancelled by a factor of M 2 in the coef- 
ficient originating from the reduction to scalar integrals. 
The purely collinear regions (lc-lc), (2c-2c) and (lc-2c) 
develop poles at several points in the parameter space of 
the rii which need to be regularized analytically and cancel 
between these three regions. 

The most complicated evaluation of the contributions 
to the Mercedes-Benz graph has to be performed for the 
(lc-lc) region with its two-fold Mellin-Barnes integral, 
especially when all propagators are present, n\ = ■ ■ ■ = 
n 6 = 1, for n 7 — 0, 1, 2. In these three cases Fbe is of order 
(M 2 /Q 2 ) -1 , as described in the previous paragraph, and 
only the (lc-lc) region contributes to the leading term. 
In addition, the integrals are finite with respect to both 
dimensional and analytic regularization, and they result in 
simply a numerical constant times (Q 2 /M 2 ). To evaluate 
these three complicated integrals, where none of the two 
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integrations can be performed explicitly due to Barnes 
lemmas (see appendix IHf . we used the following strategy 
exemplified here by Fbe(1, 1,1,1,1,1,0). After setting d — 
4 [e = 0) and applying some simplifications, equation l|32|l 
yields 



■Pi 



F BE (1,1,1,1,1,1,0) 



f dzi f d£2 
M 2 J 2iri J 2m 

— ZOO —ZOO 



r(-z 1 ) 2 r(z 1 )r(-z 2 ) 2 r(i + z 2 )r(i + Zl + z 2 ) 

1 + Z\ + z 2 



, (35) 



where the integration contours may be chosen, e.g., as 
straight lines with Re Z\ = Ylez 2 = —0.3. We performed 
the integration over z\ by closing the integration con- 
tour to the right and taking residues at the points z\ = 
0, 1, 2, . . . , to, . . ., which are given by integrals over z 2 . For 
any given m, such integrals can be evaluated with Barnes 
lemmas and their corollaries. We performed such calcula- 
tions up to order to = 100. After having understood the 
dependence of these integrals on to, we switched to "ex- 
perimental mathematics" (see e.g. [211 and EH] for ear- 
lier similar examples) and made a (successful) guess that 
the result of the integration over z 2 can be represented in 
terms of nested sums |fil| (of the argument to), in particu- 
lar sign-alternating sums. Using an ansatz as a linear com- 
bination of these nested sums, with unknown coefficients, 
we solved linear systems of equations in order to find the 
coefficients. The summation of the final series, over m, 
was quite straightforward and gave results where a value 
of the polylogarithm, Li4(i) « 0.517479, appeared: 



1 



2> ^ 4 




2 * 



■P2 



Fig. 7. Scalar graph with fermion self-energy 



£ 3 



^-7j £ 2 +(8C 3 -2 7 r 2 + |)£ 



, , M 16, 4n 16 ,„ 28 4 

128 L 4 - + — n 2 7T 2 In 2 2 tt 4 

'2/3 3 15 



3.6 Vertex correction with fermion self-energy 



(39) 



Figure |2i) shows the Feynman diagram of the vertex cor- 
rection with a self-energy insertion in one of the fermion 
lines. The reduction to scalar integrals as shown in figured 
produces the following expressions: 



F BE (1,1,1,1,1,1,0) = -^ 



7T . o _ 19 , 

— In 2 H 7r 

3 144 



F B e(1, 1,1,1, 1,1,1) =^2 
19 



-8Lu( - I - - In 4 2 



-24Li 4 ( - I - In 2 



tt 2 In 2 2 + — tt 4 - 14C 3 - tt 2 - 1 
48 



(36) 



(37) 



^(1,1,1,1,1,1,2)=-^ 



-104Li 4 (i)-H ln42 



i 13 2 , 2 9 , 247 4 31 2 73 

H 7T In 2 H tt — 63C3 7r 

3 144 SJ 6 16 



(38) 



We have checked these analytic constants by a direct nu- 
merical evaluation of the Mellin-Barnes integrals. The re- 
sult without numerator agrees with [f>2"] . 

The contributions from all relevant regions of all scalar 
integrals sum up to the vertex correction corresponding to 
the Mercedes-Benz graph in figure Efc): 



v,BE 



= \Gi 
1 

e 



1 



-CfCa 



(- 

V47T 



\M 2 



-C 2 



3£-V-H 
3 4 



F ic ( ni ,...,n 5 ) = e 2 ^ (M 2 ) 2e (Q 2 )"- 4 
d d k f d d £ 1 



fad/2 J in d/2 ( fc 2 + 2k ■ Pl )ni (k 2 + 2k ■ p 2 )™ 2 

1 



(k 2 - M 2 )™3 ((p 2 + fc + £)2)n 4 (|2 _ M 2)n 5 



(40) 



with k — k 3 , £ — k^ and n = 77,12345. For this graph, 
not every scalar product appearing in the numerator can 
be expressed linearly in terms of the only five factors in 
the denominator. But by applying standard tensor reduc- 
tion to the subgraph of the one-loop self-energy inser- 
tion (lines 4 and 5), the formally irreducible scalar prod- 
ucts may be transformed into reducible ones, so that only 
scalar integrals without numerator have to be treated. 

Due to the self-energy insertion the evaluation of the 
loop integrations is rather easy. The complete scalar inte- 
gral H4()|l with general indices rii may be expressed as an 
only two-fold Mellin-Barnes representation: 



F fc (m, . . . ,n 5 ) 



f dZ! 


f dz 2 




1 2m j 


' 2m 





e -ni7T e 2e lE p^d _ n ^ 

r(n 1 )r{n s )r{n i )r{n b ) 



r(d~n 2345 - zi) 
- 7112345 - zi) 



x r(m2345 -d + zx) r(-z 2 )r(n 3 + z 2 )r{± - m + z 2 ) 
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r(-n 3 - zi - z 2 )r{± - n 35 - zi - z 2 ) 

r(d - n 3 45 - z\- z 2 ) 
r(n 3i5 -± + Zl +z 2 ) 



r{n 234 $ 



z\ + z 2 ) 



(41) 



From this expression, the residues producing non-sup- 
pressed contributions are extracted. They correspond ex- 
actly to the contributions from the following five regions: 



(h-h) 


k r 


- Q, 


tr 




(lc-h) 


k 


pi, 


I' 




(2c-2c) 


k 


P2, 


e\ 


P2 


(h-s) 


k r 


- Q, 


lr 


- M 


(lc-s) 


k 


pi, 




- M 



r {c {rii,. 



■,n 5 ) 



These contributions are evaluated as 

Q 2 / 

r(f - n 13 )r(f - n 4 )r(f - n 5 )r{d - n 2345 ) 

r( ni )r(n 4 )r(n 5 )r(d - n 45 )r(|d - ni234 5 ) 
^("45 - ^)r(ni 234:a - d) 



r(n 2i5 - *) 



(42) 



p(lc-h) / 



,n 5 ) 



AP 

Q 1 



2— ni3+e 



k r(f - ni)T(| - n 4 )r(f - n 5 )r(| + m - n 245 ) 

r(n 1 )r(n 3 )r(n 4 )r(n 5 )r(d - n ib )r(d - n 245 ) 

x r(m 3 - f )^>45 - I) , (43) 

T 2 \ 4— n 2 345 



,(2c-2c) , \ /' Af' 

Q 



dz r(-z) 



r(n 3 )r(n 4 )r{n 5 )r(± - n x ) J 2m /> 2 - z) 

— ioo 

x r(n 24 - | - z)r(n 245 - d - z) F(n 3 + z) 

ra- ni + z)ra-n 2 + z) 

x — -, , 



(44) 




Fig. 8. Scalar non-Abelian Mercedes-Benz graph 

least one factor M 2 /Q 2 . On the other hand, the (2c-2c) 
region is of order (M 2 /Q 2 )~ 1 iin 2 = 2, n 3 = n 4 = ris = 1, 
but this inverse power of M 2 is cancelled by a factor of 
M 2 from the reduction to scalar integrals. 

For the leading order in M 2 /Q 2 , no analytic regular- 
ization is necessary. The contributions of the (h-h), (lc-h) 
and (2c-2c) regions sum up to the results for the scalar 
integrals, e.g. 



F fc (l,l,l,l,l) 



1 

2' 

+ 4C 3 



-C 2 



7T 



C 2 



F fc (l, 2, 1,1, = 



Q_ 2 

M 2 



7T 

"T ' 
i l 

~2 ~ 7 



7T 

y 



(47) 
(48) 



The whole vertex correction originating from the Feynman 
diagram in figure |2Jl) evaluates to 



Cl( 



(- 

1 

s 




C? - 3£ + ^tt 2 



v ^ 4. 171 



r(| - m 3 )r(| - n 234 )r(n 123 4 - i)r(n 5 - i) 



e -niTr g2£7B 



r(ni)r(n 5 )r(n 24 )r(d - ni 234 ) 



4 1C - S) (n 1 ,...,n 5 )=(^- 



M 7 



2 \ 4-ni 35 



g— ni7r g2e7E 



(45) 



x i>r - » 24 )r(f - m)r(ni 3 - f )/> 5 - f ) 

r(ni)r(n3)r(nB).r(g-n34) 

The contribution of the (h-h) region is known from the 
massless diagram 53 54 ,55,64 . In the reduction to scalar 
integrals only parameters raj with n 2 < 2 and m < 1, 
i = 1,3,4,5, are involved. Therefore the contributions of 
the (h-s) and (lc-s) regions are always suppressed by at 



3.7 Vertex correction with non-Abelian 
Mercedes-Benz graph 

The Mercedes-Benz graph in figure 01) is of pure non- 
Abelian nature due to its three-gauge-boson vertex. The 
corresponding scalar integrals are illustrated in figure |H1 
and defined as follows: 

F B EC,(ni, .. . ,n 7 ) = e 2 ^- (M 2 ) 2e (Q 2 )"-"- 4 
d d k f d d £ (2p 2 • fc)™ 7 



in d/2 J in d/2 ( £ 2 _ 2pi . £)ni ((2 _ 2 p 2 . £)r 
1 

(k 2 - 2 Pl ■ fc)™3 (p - M 2 ) n * (k 2 - M 2 ) n * 

1 



((jfe - I) 2 - Af 2 )™« ' 



(50) 
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with k — fcg, I — ki and n = "123456- This definition 
is the same as for the Abelian Mercedes-Benz graph in 
figure Eland equation (|2"?|l . except for the distribution of 
the masses in the propagators. Also the list of relevant 
regions is similar, but now there is the (s'-h) region instead 
of the (us-2c) region: 



(h-h) 


k r 


-Q, 


£r 




(lc-h) 


k\ 


pi, 


I' 




(h-2c) 


k r 


- Q, 


I 


P2 


(s'-h) 


k e 


~M, 






(lc-lc) 


k\ 


Pi j 


t 


pi 


(2c-2c) 


k\ 


P2, 


i 


P2 


(lc-2c) 


k\ 


Pi j 




Pi 



The (s'-h) region is of order (M 2 /Q 2 ) 2 ~™ 6+£ and there- 
fore suppressed with respect to the (h-h) region, as 
rii < 1 (i = 1,...,6). The (2c-2c) region is of order 
(M 2 /<3 2 ) 4 "™ 2456+ ™ 7 and suppressed for n 7 > 0; it is only 
evaluated for n 7 — 0. 

The leading contributions of regions with k 3 ,k§,kQ ~ 
Q are identical to the corresponding contributions of the 
Abelian Mercedes-Benz graph: 



p( h " h ) (r, 



BEG, ( 
p (h-2c). 



,n 7 ) 



^BE l n l' 



,n 7 ) , 

BE( n 7 ) — Fg E \m, . . . , n 7 ) . 
The contributions of the other regions are given by 

2 \ 2-ri 35 +e 



(51) 
(52) 



p (lc-h). 



.,m) 



Q 2 



g— nlTT g2E7B 



r(f-n 3 )r(f - n 146 )/>35 - f ) 
r(ni)r(n 2 )r(n 3 )r(n 5 )r(n e )r(d - ni246) 



dz r(-z)F(| - n 24 6 - z)-r(n 6 + z)r(n 37 + z) 



2m 

X r(ni246 
^(lc-lc) / 

il23<™7 
\ " 



n 7 



g + z). 



(53) 



M 2 

Q~ 2 



n 7 \ 



r(n l7 - ii 2 ) 



ii ifi^ >o * l! * 2 ' * 3! ~ il23 ^' r ( n i) r ( n 3) r ( n i)r(n 5 ) 

zoo ZOO 

dzi f dz 2 , . , . 



r(n 6 )r(f - n 2 + n r ) 

— zoo —zoo 

x T(n 5 + zi)r(| - n 3 + n 7 - i 123 + Z\)r{m + h + z 2 ) 
r(n li67 - i 123 - I + z 2 )r(ni34 56 - d + z 2 ) 



-T(ni34 



r{n li7 -i 2 + z 2 ) 

*13 - I - 21 + 2 2 ) 



-T("14567 - «123 - 2 + z l + z 2) 
"-14 + i 2 + Zi - Z 2 ) 



7^(2 - "4 + n 7 — «1 + Z\ — z 2 ) 
x r(| — n 2 4 + n 7 — ii + zi — z 2 ) , 



(54) 



,-,(2c-2c)/ . 

f bec a {ni,...,n 6 ,n 7 = 0) = . 

7"(n 2 - n i3 ) 



M 2 



2 \ 4-n 2 456 



g -niw g2£7-E 



r(n 2 )r(n4)r(n 5 )r(n 6 )r(|-n 13 ) 7 27ri 

— zoo 

r{n 2A - f - z)T(n 5 + z)r(n 6 - n 3 + z) 

± 

r(n 56 -n 3 + 2z) 
x r(| - n 2 + z)J> 56 -i+z) 



■n-z) 



> 2 

-,(lc-2c). 



2 / ' 

2 \ 4 — H2345 



(55) 



g-nwr g2e7B 



T(n 2 - ni6)r(n 37 - n 6 )r(| - n 2 )r(f - n 3 ) 

x r(n 2 )r( Il3 )r( )1 4)r( ) i5)r(f - n M )r(f - n 6 + n 7 ) 

x T(n 24 - f )r(n 35 - f ) . (56) 

The evaluation of this non- Abelian vertex graph is more 
complicated than in the Abelian case, mainly due to 
the appearance of three massive propagators. The com- 
plete summation of the infinite number of residues in the 
Mellin-Barnes integrals (|54|l and Q55JI is quite intricate. 
As the calculation of the four-fermion amplitude demands 
the result of the form factor only to N 3 LL accuracy, we 
have refrained from calculating the non-logarithmic con- 
stant in the non- Abelian corrections (cf. the beginning of 
section^ . Therefore we have only extracted all logarithms 
ln(Q 2 /M 2 ) from the integrals. 

The (lc-h) region has been evaluated in the usual way 
as described in the previous sections. From the (c-c) re- 
gions, i.e. (lc-lc), (2c- 2c) and (lc-2c), the logarithmic con- 
tributions have been isolated. The expressions for the re- 
gions depend on M 2 /Q 2 only through a prefactor of the 
form (M 2 /Q 2 ) m+X , where m is an integer and x is made 
up of regularization parameters (like s) tending to zero. 
Logarithms ln(Q 2 /M 2 ) arise only when poles in the regu- 
larization parameters appear, e.g. 



M 2 Y 1 1 . (Q 2 , , , 
Q 2 -) x^- X - ln {W 2 ^ +0{x) 



As the exponents of these prefactors in the contributions 
1154( 1 156|l of the (c-c) regions involve only the parame- 
ters rii and not e, only poles originating from the analytic 
regularization may give rise to logarithms. A thorough 
analysis of the Mellin-Barnes integrals shows that such 
poles only appear in the following seven integrals: 

^(-1,1,1,1,1,1,0), 

f£-£(o, 1, 1, 1, 1, 1, m) with n 7 = 0, 1, 2, 

f£-£ (1,1, 0,1, 1,1,0) and 

ifcJC 1 : 1. 1, 1. 1. 0, n 7 ) with n 7 = 0, 1. 

When closing the integration contours in the Mellin- 
Barnes integrals, it is sufficient to take those residues 
which are responsable for the poles. In most cases these 
are only a finite number of residues. Only the integrals 
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F^ A (0,1, 1,1, 1,1,0) and (1, 1, 0, 1, 1, 1, 0) need 

the summation of an infinite number of residues. Such 
summations can be transformed to infinite series like 



1 



^ ( 2m ) \ >! - _>//> 

m=0 V m , 

1 



1 



1 



l + 2m 



1 



3 + 2m 1 + 2m 



+ Si(2m)-5i(m) , (57) 



where Sx(m) — J2iLi \ 1S a harmonic sum. Recently there 
has been a lot of progess in solving summations like l(5T]) 
analytically (see e.g. IrJSIIfjcJ] '). But still not all possible 
cases are covered, and the transformation of a given ex- 
pression into a series where the solution is known can be 
quite cumbersome. 

On the other hand, the series in (|57|l is converging very 
fast. The numerical summation of the first 300 terms ap- 
proximates the series with an accuracy of more than 100 
decimal digits. This enables us to use the following method 
(see e.g. jSZ])- An ansatz is chosen as a linear combination 
of analytical constants like 7r 2 , £3, In 4 2 etc. with unknown 
rational coefficients. The determination of the coefficients 
starting from the numerical result is performed by the 
PSLQ algorithm 68 69 70 . We have used an implemen- 
tation |71j of PSLQ in Fortran with multiprecision arith- 
metic The series above in ifBTt has hereby been 
identified with the analytical expression 4v / 3Cl2( : | 
where C^f) ~ 1.014942 is a value of the Clausen func- 
tion. 

In addition to the logarithms, we have calculated in 
a purely analytical way the complete set of poles in e in 
order to control the cancellation of ultraviolet and infrared 
singularities. 

The contributions from all regions sum up to the re- 
sults of the scalar integrals, from which we quote the two 
most complicated ones: 

FbkcAO, 1, 1, 1, 1, 1, 0) = -^£ 4 - y£ 2 + ^( 3 £ , (58) 

Fbec, (1, 1, 0, 1, 1, 1, 0) = J (-^C 2 ~ y) + \& - C? 

+ 2V3C1 2 (|) £, ( 59 ) 

where non-logarithmic terms of order e° have been omit- 
ted. The integrals Fbec a (1j n 2, 1, 1, 1, 1, n 7) are of order 
(Af 2 /Q 2 ) -1 , with the only contribution coming from the 
(lc-lc) region, and the inverse power of M 2 is again can- 
celled by a factor of M 2 from the reduction to scalar in- 
tegrals. These integrals, however, produce neither loga- 
rithms C = ln(Q 2 1 M 2 ) nor poles in e and do therefore 
not contribute to the result in N 3 LL accuracy. The result 
for the vertex correction corresponding to the non-Abelian 
Mercedes-Benz graph in figure |3Ji) is as follows: 



/ a \ 2 / /1 
bvc a = C f C a {—) l w 



1 





Fig. 9. Scalar graph for non-Abelian vertex corrections with 
loop insertions 



2 I 6 2 1 



W3Cl,( ?H C3-^ + ?U 



+ O(e"£ ) +0(e) + O 



[Ml 

\Q 2 



(60) 



g 3.8 Non-Abelian vertex corrections with loop insertions 



This section treats all vertex diagrams from figure|3where 
a self-energy loop has been inserted in the gauge boson 
propagator: the gauge boson loop in figure |2t>) , the ghost 
field loop in figure |3J:) and loops involving the Higgs and 
Goldstone bosons in figures |3Ji) , 03) and |3f ) . Care must 
be taken in the interpretation of the Feynman rules (ap- 
pendix not to forget the factor (—1) for the loop of 
the anticommuting ghost fields and the symmetry factor 
1/2 for the loops with two gauge bosons or two Goldstone 
bosons. 

Additional contributions from "tadpoles" , where a 
loop of only one gauge, Higgs or Goldstone boson is at- 
tached to the gauge boson propagator via a vertex with 
four fields, are omitted here because they are cancelled ex- 
actly by the corresponding contributions from the renor- 
malization of the gauge boson mass (see section 0J. 

For the Higgs boson mass Mjj we use the approxima- 
tion Mji = M\y, which facilitates the loop calculations. 
The form factor depends on the Higgs mass only in N 3 LL 
accuracy, i.e. via the coefficient of the linear logarithm, 
and the higher powers of the electroweak logarithm are not 
affected by changes in the Higgs mass. We have checked 
explicitly by evaluating the Higgs contributions for the hy- 
pothetical case Mh = (see the discussion of the results 
in section |SJ) that effects due to a wrong value of the Higgs 
mass are indeed negligible. 

As in our approximation (and using the Feynman- 
't Hooft gauge) all particles running in the self-energy 
loop have the same mass M — M\y, the vertex correc- 
tions of this section share the same set of scalar integrals, 
which are illustrated in figure El and defined in the follow- 
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ing equation: 

F Wc (n 1 , ...,n 6 ) = < 
d d k f d d 



2e lE ( M 2)2e (Q2)n-4 
1 



ilT d/2 J in d/2 ^2 + 2 Pl . (A: 2 + 2p 2 ■ fc)" 2 
1 

X (fc 2 - M 2 )"3 (£2 _ M 2)n 4 ((£ + £)2 _ M 2)n 5 (£2)^ ' 

(61) 

with k = k 3 and I = k 4 . The additional massless prop- 
agator corresponding to the parameter «6 is introduced 
when performing a tensor reduction on the self-energy 
loop (lines 4 and 5) in order to eliminate the scalar prod- 
ucts in the numerator of the integral. The presence of both 
a massive and a massless propagator with the same mo- 
mentum &3 = fcg — k could, of course, be avoided by par- 
tial fractioning. But this would produce factors of 1/M 2 , 
complicating the expansion in M 2 /Q 2 . So we remained 
with both propagators in the scalar integrals l|bT|) . In or- 
der to avoid ambiguities in the reduction to scalar inte- 
grals, we fixed n 3 = 2 and cancelled factors of k 2 in the 
numerator exclusively with the sixth propagator. 
In general the following regions are relevant: 



(h-h) 


k r 


-Q, 


Ir- 




(h-s) 


k r 


-Q, 


tr- 


- M 


(h-s') 


k < 


-Q, 


h 


~ M 


(lc-lc) 


k 


Pi, 


*\ 


Pi 


(2c-2c) 


k 


P2, 


l\ 


P2 



But since the (h-s) regions is of order (M 2 /Q 2 ) 2 ~ ni+£ and 
the (h-s') regions is of order (M 2 /Q 2 ) 2 ~ ns+E , they are 
both suppressed with respect to the (h-h) region for all 
relevant cases. The contributions from the other regions 
can be expressed as follows: 



F 



(h-h) 



2? 



g-niir e 2e lE 



M 2 ' 

wc '{ni,...,n 6 ) = \ ~q2 

T(f - n 4 )r(f - n 5 )r(d - n 13456 )r{d - n 23456 ) 
r(n 1 )r(n 2 )r(n4)r(n 5 )r(d - rc 45 )F(§d - ni2345e) 



x r(n 45 - |)r(ni23456 - d) 



(62) 



F 



(lc-lc) 



Wc 



("1 



,n 6 ) 



M 2 



2 \ 4 — 7113456 



Q 2 

r(ni - n 2 ) 



g -row g2e7£ 



- : n-z) 



, dz 

r( ni )r(n 3 )r( ni )r(n 5 )r^ - n 2 ) J 2™ 

— ioo 

F(n 136 -|- z)r(n 4 + z)r{n 5 + z)r{± - n 16 + z) 



r{n 45 + 2z) 



x r(n 45 ~ i+z). 



(63) 



originating either from the reduction to scalar integrals 
or from the Feynman rules, when two WWH -vertices are 
present. 

As for the non-Abelian Mercedes-Benz graph in the 
previous section, we have only extracted the logarithms 
and the poles in e. The (c-c) regions (lc-lc) and (2c-2c) 
produce logarithmic terms for n± = n 2 = 1. For most of 
the (c-c) contributions, the evaluation of a finite number 
of residues in the complex z-plane is sufficient. Only the 
two cases F$f(l, 1, 2, 1, 1, 0) and F^ c c) (l, 1, 2, 1, 1, -1) 
demand the summation of an infinite number of residues. 
We have solved these two summations numerically and 
found the corresponding analytic expressions with the help 
of the PSLQ algorithm. In all cases the extraction of the 
poles in e required only a finite number of residues. 

We quote the results for the two most complicated 
scalar integrals: 



F Wc (l, 1,2, 1,1,0) = -^ 



F Wc (l, 1, 2, 1, 1,-1) = X - (- l -C 2 



1 



C 



2 



2 \C. 



(64) 

3 

(65) 



The vertex corrections corresponding to the Feynman di- 
agrams in figures|3|3)-|3f ) are obtained by inserting the re- 
sults for the scalar integrals into the expressions returned 
from the reduction of each diagram. 

The vertex corrections with the non-Abelian gauge bo- 
son and ghost field loops, figures Ot) and|2t), have been 
evaluated together. Their sum is 



F 



v,WWcc 



= CfCa 

1 

e 



( a 
V4^ 



\M 2 



2 c 



Si 



5 r2 49 10 2 337 

Y + Y c -Y n ~l2 



76 



-4V3C1; 



859 
~L8~ 



10 3 
9 



C 



O{e £°) + O(e) + o( 



M 2 

Ko 1 



(66) 



The vertex correction in figure |3JI) with gauge and Higgs 
boson in the loop insertion contains two factors of M from 
the two W WH- vertices, but these are cancelled by factors 
1/M 2 in the results of some of the scalar integrals. So this 
vertex correction is not suppressed: 



,WH 



for\ 2 ( 



V47T 



M 2 











I' 





For symmetry reasons, F^" 20 -* can be obtained from 



Fwc" 1C ' by exchanging n\ «-> n 2 . Some of the (lc-lc) and 
(2c-2c) contributions are of order (M 2 /Q 2 )^ 1 , but this in- 
verse power of M 2 is always cancelled by a factor of M 2 



2V3C1 2 



3 C 



O(e C°) + 0(e) + O 



[Ml 



(67) 
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The two vertex corrections in figures Ot) and U) with 
Higgs and Goldstone bosons in the loop insertion yield 
the same result (for Mh = M): 



16 16 24 64 



-— C? 
24 



17 



■O(e C°) + O(e)+of 



3 



19 
96 



C 



(68) 



The contributions involving Higgs and Goldstone bosons 
are only valid for the spontaneously broken SU(2) model, 
they cannot be transformed e.g. to a U(l) model simply by 
setting other values for Cf, Ca and Tp. We have therefore 
written these contributions with the Casimir operators al- 
ready replaced by their SU{2) values. 



4 Renormalization contributions 

Section has treated the evaluation of the vertex correc- 
tions which contribute to the Abelian vector form factor. 
These have been performed with Feynman rules originat- 
ing from the unrenormalized Lagrangian. Therefore the 
contributions due to the renormalization of the fields (sec- 
tion l4. lfl , the coupling constant fsection l4~2")l and the gauge 
boson mass (section |4.3|) have to be added. 



4.1 Field renormalization 

The renormalization of the two fermion fields in the 
Abelian vector current requires the multiplication of the 
vertex corrections by a factor of Zf, where {Zf) 1 / 2 is the 
fermion field renormalization constant. On the other hand, 
Zf is determined by the fermion self-energy corrections S 
at on-shell momentum p 2 = (for massless fermions) . In a 
perturbative expansion, the field renormalization constant 
is Zf = 1 + Si + £2 + 0(a 3 ) and the vertex corrections 
are F v — 1 + F v _i +-F V ,2 + 0(a 3 ), where the indices 1 and 2 
indicate the one- and two-loop contributions, respectively. 

The total Abelian vector form factor up to order a 2 
can be written as 

F = F v ■ Z f 

= 1 + F Y ,i + £1 + F Vi2 + E 2 + F Vtl Si + <D(a 3 ) . 

O(a) 0(a 2 ) 

The two-loop contribution to the form factor is therefore 
given by 

F 2 = F v ,2 + £2 + F vA £i , (69) 




Fig. 10. Fermionic self-energy correction 




(a) 



(b) 



(c) 

Fig. 11. Abelian self-energy corrections 



where -F v ,2 is made up of the contributions calculated in 
section |3| The one-loop corrections are well known: 




C + 2( 3 + -ir 2 



-c 4 + -c 3 

12 2 



7i 

12 



4 \£ 2 



4 -- +16 J £ -ilr 4+ T^ 2 



28 



(70) 



2?i = £i(p 2 = 0) = C F 



0(e 3 ). 



,2 \ £ 



4tt V M 2 



1 1 

7+2 



" IT 2 1" 


to 


12 + 4_ 


1 



1 A 7T 2 1 

3 C3+ 24 + 8 



(71) 



The terms proportional to e and e 2 are needed when F Vi i 
and S\ are multiplied by other one-loop contributions 
containing 1 /e-poles from ultraviolet singularities or 1/ e 2 - 
poles from mass singularities. The sum F\ — F Y) i+Si con- 
stitutes the one-loop form factor, which is finite at e = 0. 

The two-loop self-energy corrections S2 originate from 
the Feynman diagrams in figures I1UI 1771 and ^| As for 
the vertex corrections, "tadpole" diagrams are omitted 
because their contributions are cancelled by the renormal- 
ization of the gauge boson mass. 

The self-energy amplitudes S are quadratic matrices 
both in the spinor and in the isospin space. For massless 
fermions of momentum p they are of the form 



S(p) = -ifflS( P 2 



(72) 
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(a) 



(b) 





(c) 



(d) 




(e) (f) 
Fig. 12. Non-Abelian self-energy corrections 



where 1 is the unity matrix in the isospin space. The self- 
energy correction E may be extracted from the ampli- 
tude E by the projection 



E 



4Np 2 



Tr 



(73) 



where N — 2 for 5(7(2) and the trace runs over the spinor 
and the isospin indices. The projection requires p 2 ^ 0, 
whereas we need the self-energies at p 2 = 0. We have 
therefore calculated the loop integrals for an infmitesi- 
mally small, but finite p 2 . By performing the limit p 2 — ► 
before any expansion in e, no logarithms ln(p 2 ) appear 
and the first two coefficients of a simple Taylor expansion 
of the integrals with respect to p 2 are sufficient. The con- 
tribution of every Feynman diagram to the trace in l|73|l 
is proportional to p 2 , so no inverse power l/p 2 is left. 

The reduction of the self-energy diagrams to scalar in- 
tegrals using (|757 and the calculation of the integrals was 
performed similarly to the vertex corrections, fn fact, the 
evaluation of the self-energy corrections is much simpler 
because they do not depend on Q 2 , only on M 2 , and no 
expansion by regions is needed. We do not quote further 
details of this calculation and list only the total results of 
each Feynman diagram. 

The fcrmionic self-energy correction in figure has 
already been calculated in |2*2*] : 



C F T F n f (^Y(^ 



V47T 



M 2 



2 s 



0{e). 
(74) 



The Abclian contributions to the self-energy correction 
originate from figure ITTk ). 



E 



Tl 



cl 



1 



-CfCa 



from figure HTbl. 



f- 

3 

2e ' 



2 s 



7T 
T 



M 2 
7 

+ 4 



2? 



0(e), (75) 



E< 



T2 



\M 2 J 



1 

2^ 



i 



7T 

12 



0(e). 



(76) 



and from figure mb ) . which yields just the square of the 
one-loop correction l|7T|) . Only the C|.-part of lj75|) belongs 
to the Abelian corrections, the (TpC^-part contributes to 
the non- Abelian corrections. 

The other non- Abelian contributions have been evalu- 
ated from figure H2k,). 

E t1Ca = C f C a (^-) (jL^ 5 £ 2 (-A~|) 

+ O(e ), (77) 

from the sum of the diagrams with gauge boson and ghost 
field loops, figures Eb) and 112b). 



EwWcc = CfCa y^r— 



(V 

\M 2 



S 2 £ -^- + (D(e a ), (78) 



from figure I12t with gauge and Higgs boson in the loop 
insertion, 



Ewh — I — 

47T 



4s 



(79) 



and, with identical results, from figures!*^ - ) andll2t' % ) with 
Higgs and Goldstone bosons in the loop insertion, 



E 



H4> 



( a 
V47 



JJ_ 

M 2 



2 s 



S 2 e ~ £ +O(e ). (80) 



As for the vertex corrections, the Higgs and Goldstone bo- 
son contributions have been calculated with the approx- 
imation Mh = M and are only valid in a spontaneously 
broken SU(2) model. The evaluation of all non-Abelian 
contributions has been limited to the poles in e because 
the non-logarithmic finite term of order e° has already 
been neglected in the calculation of the corresponding ver- 
tex corrections. 



4.2 Coupling constant renormalization 

According to the prescription of the MS scheme, the un- 
renormalized coupling constant abate is replaced by the 
renormalized coupling a via 



(81) 
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where (3q is the one-loop coefficient of the renormaliza- 
tion group /3-function. /3q gets a non-Abelian contribution 
proportional to Ca, a fermionic contribution proportional 
to n/ and a Higgs contribution [741175] : 



(82) 



As mentioned above, the loop calculations have been per- 
formed using the unrenormalized Feynman rules. Intro- 
ducing now the renormalized coupling constant and mass 
instead of the bare quantities does not change the two- 
loop results at order a 2 . But the coupling and mass in the 
one-loop result have to be regarded as the bare parameters 
and must be replaced by the renormalized ones. 

By applying the substitution (|81|l to the one-loop form 
factor from equations (|70l) and (|71|l . we get additional con- 
tributions of order a 2 , namely 



11 



,2 n 22 2 77 

C? - 1LC+ — 7T 2 1 



9 



11 



11 



-£ 2 + 



11 



„ 22 , 



11 



517 



— ) + 0{e)+0[-^ 



12 



A7 2 



(83) 



2 \ e 



m 2 y 



4., „ 8 „ 14 

--£ 2 +4£- -7T 2 - y 

4 2 32\ 8, 
-9* + ¥ £+ 3 Ca " 



-£ J - 2£ 2 



3* 



Higgs ^ 



M 2 



\Q 2 

2 \ e 



(84) 



-i£ 2 + ^-^-l 
8 8 12 16 



i£3 

24 



1 



-^ £2+ (-i4 +1 J £+ 4 C3+ 16 




(85) 



4.3 Mass renormalization 

The relation between the bare gauge boson mass Mbare 
and the renormalized mass M is determined by the gauge 
boson self-energy corrections, which have the form 

]J^ ab (k) = iS ab g^k 2 n(k 2 ) + terms oc k^k" (86) 



at momentum k. In the on-shell scheme, the square of the 
physical, renormalized mass is defined to be the real part 
of the pole of the propagator. At one-loop, the relation 
between Mbare and M becomes 



M 2 arc = M 2 [1 - Re7Ti(M 2 )] + 0(a 2 ) 



(87) 



where 77i is the one-loop contribution to 77. This relation 
leads to the following substitutions in the one-loop result 
(HI) and (El): 



Af 2 



jfiy 

A7 2 I 



[1 + e Re77!(A7 2 )] + 0(a 2 



C n -» C n + nC l ~ x Re77 1 (A7 2 ) + 0(a 2 ) . 



(89) 



with C = ln(Q 2 /A7 2 ), producing additional contributions 
of order a 2 . 

The one-loop gauge boson self-energy receives contri- 
butions from a fermion loop, 



U <M*\ T 01 ( a ( 4 20 4 ' 

ll n AM = ifii/ — — k j e w 

A ' 3 4tt VM 2 / e V 3e 9 3 

+ 0(e), (90) 

from the non-Abelian gauge boson and ghost field loops, 



77 Wcc (A7 2 ) = CU £ (j^ 



17 4tt 82 



+ (91) 



from the loop with gauge and Higgs boson 

" .2 X * / ! 



+ 0(e), (92) 
and from the loops with Higgs and Goldstone bosons, 



n m {M 2 ) = n^(M 
— (^-\ s 

4tt I A7 2 1 



5 

12e 



4\/3 



17 
18 



0(e). (93) 



The self-energy diagrams with "tadpoles" have been omit- 
ted. They do not depend on the momentum of the gauge 
boson, so their contribution to the mass renormalization 
cancels exactly the corresponding vertex correction and 
field renormalization diagrams which have already been 
dropped out before. 

Applying the substitutions l|88|l and (|89fl to the one- 
loop form factor, the self-energy corrections (|9()|) - (|93[) pro- 
duce the following contributions to the two-loop form fac- 
tor: 



AF% = C F T F n f 



40 „ 




Bernd Jantzen, Vladimir A. Smirnov: The two-loop vector form factor in the Sudakov limit 



17 



AR 



M 

WWcc 



A T? M 



34 „ 
- y £ + 17 



8tt 164 \ „ n 

-=-—)£- 4V37T - — TT 2 

\/3 9 J 3 
/ M 2 



-C-- 

2 4 




(95) 
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5 Results and discussion 
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(99) 



For the non-Abelian contributions proportional to CfCa 
the remaining part of Fv.np, -Fv.be and I?ti is considered 
together with the purely non-Abelian results from equa- 
tions lEU, (MJ), (El, (|H3I» and JHHl): 



2,Cf Ca 



NP 



2 fv.BE + -£ti 



*Fg A 



A J? M 



233 
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4 ^ cl2 (3) + ^-y C3 + y^ 



>/3 

i??) £}+O(£ ) 



(100) 



The individual results have been presented in the previous 
sections so that we can now add them together. Accord- 
ing to the MS prescription, the factor S s — (47r) e e~ e7E 
is absorbed into /x 2£ by a redefinition of fi, and we have 
chosen /i = M so that the whole prefactor (/i 2 /M 2 ) e S e 
or (/i 2 1 'M 2 ) 2e S 2 is replaced by 1. The dependence of the 
form factor on /1 can easily be restored by looking at the 
running of the coupling a, parametrized by /3o, in the one- 
loop form factor. We give the results in d = 4 dimensions 
(e = 0) in the Sudakov limit Q 2 > M 2 . 

The fermionic contribution to the Abelian vector form 
factor is obtained from equations (0, (|74|l . 184(1 and 
lHU) [SI: 



Fl,n s — Fv,n f + £n f + AF" 



16 2 

+ 7T 

27 



Utt 



-C 



AF M 
38 



34 



115 
~~9~ 



(98) 



with C — ln(Q 2 /M 2 ). For the Abelian contributions only 
the Cp part of F Vj np, -Fv,be and Xti is considered. The 
vertex corrections Fv,be and -F v ,fc have to be counted twice 
because two horizontally mirrored diagrams exist for each 
of these. The result follows from equations (|14[l . I|26|) . (|39|) . 
P9"|). J70J), CD, G3> and E7j: 



F 



2,Cf 



-Fv,la + -Fv,np|c|, + 2-F v ,be|c| + 2-FV.fc 

+ -^Tllcf. + -^T2 + (^l) 2 + F v ,l£l 



The Higgs contribution results from equations (|67|l . I|68(l . 
(El, ®, (J2U, © and JHIIl: 

i*2, Higgs = -Fv,VKH + F v .H<t> + Fy^tj, + Sv/H + ^Hct> + ^00 

+ AF& ggs + AF^ H + AF™ + 

1 r- TT 2 23 

+ 4^+ 24 + 16 



£|+O(£ ). (101) 



The two non-Abelian contributions F2.c F c A and i 7 ^, Higgs 
depend on the Feynman-'t Hooft gauge in which they have 
been calculated. Only their sum is gauge invariant: 

•F2,CVCA+Higgs = F 2 .c F C A + ^i.Higgs 

/M 2 (43 907\ 
\4ttJ {24 \ 2 48 J 



f ^ - 44(3 + |» 2 



(102) 



where the values Cf = 3/4 and Ca = 2 for the SU(2) 
gauge group have been used. Adding all contributions to- 
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gether, the two-loop form factor is given by |33| 

F 2 = F 2 ,n s + F 2 ,c 2 F + F 2: c F c A +m se s 




+ 0(£°). (103) 

The coefficients of the first three logarithms £ 4 , £ 3 and £ 2 
agree with the NNLL prediction of the evolution equation 
approach jl lU12j . The coefficient of the linear logarithm is 
a new result. 

Let us have a look at the numerical size of the coef- 
ficients in the individual contributions. For the fermionic 
contribution we set rif = 6 for 3 lepton and 3x3 quark 
doublets from which only the left-handed degrees of free- 
dom couple to the gauge bosons. 

F 2 , n , « (£^) 2 ( -1.0£ 3 + 9.5£ 2 -26£ + 42) , 
F 2 C 2 p sa (£-) 2 (0.3£ 4 -1.7£ 3 + 8.2£ 2 -ll£ + 15), 

i^.CjrCji+Higgs 

«(^) 2 ( 1.8£ 3 -14£ 2 +43£ + ...) 

(104) 

We notice that all three contributions show a similar pat- 
tern of coefficients with alternating signs and growing size. 
At a typical energy in the TeV range, Q = 1 TeV, using 
M = 80GeV and a/(4?r) = 0.003 as rough values for the 
weak interaction, the individual logarithmic terms have 
the following numerical size in per mil (1/1000): 

£ 4 £ 3 £ 2 O £° 
F 2 , n/ -> -1.2 + 2.2-1.2 + 0.4, 

F 2 ,c 2 F ->■ + 1.6 - 2.0 + 1.9 - 0.5 + 0.1 , 
F 2 ,c F c A +m sss -> +2.1-3.2 + 2.0 + ... (105) 

The pattern of growing coefficients with alternating signs 
produces large cancellations between the terms of differ- 
ent powers of logarithms and also between F 2 _ nf , F 2 c i 
and F 2i c F c A +Higgs- In each line of (|105f) . the largest term 
is reached at the quadratic or (for F 2 C ^) already at the 
cubic logarithm. The linear-logarithmic term is less signif- 
icant and, at least for the fermionic and the Abelian part, 
the non-logarithmic constant is again smaller by a factor 
of 3 or more. For the sum of the three contributions, 

£ 4 £ 3 C 2 C 1 
fW + 1.6 -1.0 + 0.9 + 0.3, (106) 

the logarithmic terms are monotonically decreasing in size 
already from £ 4 on. Due to the cancellations between the 
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Fig. 13. Two-loop contribution to the Abelian vector form 
factor F2 in successive logarithmic approximations, using the 
values M = 80 GeV and q/(4tt) = 0.003 

individual contributions in I|1U5|) , the non-logarithmic con- 
stant of Fj^. is larger than the total linear-logarithmic 
term in (|106f) . But the logarithmic terms in all contribu- 
tions and in the total form factor are getting significantly 
smaller from the linear logarithm on. So we do not expect 
the neglected non-logarithmic constant of the total result 
to be larger than the total linear-logarithmic term. This 
leads us to the conclusion that the N 3 LL result with all 
logarithmic terms approximates well the full result. 

Figure ^| illustrates the behaviour of the successive 
logarithmic approximations, starting from the LL approx- 
imation with only the £ 4 term and adding one after the 
other the smaller powers of logarithms. 

The result presented here relies on the approximation 
that the Higgs mass is equal to the gauge boson mass, 
Mh = M. In order to investigate the dependence of the 
form factor on the Higgs mass, we have also calculated the 
Higgs contributions in the hypothetical case of a vanishing 
Higgs mass, Mb = 0. Then equation IjlUll) becomes 




(107) 

Only the coefficient of the linear logarithm differs between 
equations (|101fl and (|107fl . The coefficients of the cubic 
and quadratic logarithms are the same, they do not de- 
pend on the Higgs boson mass and have already been de- 
termined in the evolution equation approach |1 lU12j . By 
setting Mjj = 0, the coefficient of the linear logarithm of 
-p2,cvCA+Hi g gs in (|104J) numerically changes from 43 to 45, 
and the contribution of this term in (|105|) is shifted from 
2.0 to 2.1. So the variation of the N 3 LL form factor be- 
tween the two cases Mh = M and Mh = is smaller than 
the total linear-logarithmic contribution by a factor of 3. 
On this basis we expect the deviation of the form factor 
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with the true Higgs mass from our result to be comparable 
to the neglected non-logarithmic constant. 

Altogether we estimate the accuracy of our form factor 
result to be of the order of the linear-logarithmic contri- 
bution, i.e. about half a per mil with respect to the Born 
result. 

The result for the Abelian vector form factor presented 
in l| 103(1 has been combined in with the reduced 

amplitude from equation QJ in order to obtain the four- 
fcrmion scattering amplitude in the spontaneously broken 
SU(2) model in N 3 LL accuracy. In addition predictions 
for the electroweak model have been obtained by separat- 
ing the infrared-divergent electromagnetic contributions 
(cf. appendix [Dj) and by expanding in the mass difference 
between the W and Z bosons. For a discussion of this pro- 
cedure and of the accuracy of the electroweak corrections 
we refer to [351131] . 



6 Summary 



Vertices involving four fields and vertices without a 
gauge boson do not appear in our present calculation and 
are omitted here. 



Gauge boson coupling to fermions 




Gauge boson self-coupling 




a be 



9 Pf J.( k 3 



k 2 ) P 



In the present paper we have discussed in detail the cal- 
culation of the two-loop corrections to the Abelian vector 
form factor in a spontaneously broken SU (2) model. The 
result was obtained in N 3 LL accuracy and contains all log- 
arithmically enhanced terms. It enables the derivation of 
electroweak corrections to four-fermion processes with an 
error of a few per mil to one percent, thus coping with the 
expected experimental accuracy at a future linear collider. 
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Gauge boson coupling to ghost fields 



- c 



Gauge boson coupling to Higgs and Goldstone bosons 



w 



k H 



H =igM gil „5 



H 



ah 



k <t>)^ 



A Feynman rules 

This appendix lists the Feynman rules of the vertices 
which are needed for the calculation of the form factor, 
as they follow from the Lagrangian of the spontaneously 
broken SU{2) gauge model described in section 

The gauge boson fields of mass M = M\y are W£, 
a = 1, 2, 3 (with Lorentz vector index /i). To each W a cor- 
responds a ghost field c a (and antighost c a ) and a Gold- 
stone boson 4> a , one of the unphysical components of the 
Higgs doublet. In the Feynman-'t Hooft gauge used by us, 
there is M c = = M\y- The physical Higgs boson H 
has the mass Mjj- Finally, ip denotes a fermion (lepton or 
quark) doublet of Dirac spinors, and g is the weak SU(2) 
coupling. 



K2 



In contrast to the other vertices above, which can be used 
in any SU (N) gauge model, the couplings involving Higgs 
and Goldstone bosons are only valid for the spontaneously 
broken SU(2) model. 



B Expansion by regions 

The asymptotic expansion of Feynman integrals in limits 
typical of Euclidean space is given by well-known prescrip- 
tions as a sum over certain subgraphs 



20 



Bernd Jantzen, Vladimir A. Smirnov: The two-loop vector form factor in the Sudakov limit 



However, the Sudakov limit we are dealing with is typical 
of Minkowski space. Still for some special cases, similar 
graph-theoretical prescriptions were obtained |81II82II57| . 
In particular, as it was shown in |57| . they can be applied 
to expand the planar two-loop diagram of figure^) in the 
Sudakov limit. 

The bad news is that for other relevant diagrams, 
such as the non-planar and the Mercedes-Benz diagrams, 
graph-theoretical prescriptions are not available. The good 
news is that one can apply here (and for any other limit) 
the strategy of expansion by regions [H3, 58 ,84 52 which 
consists of the following prescriptions: 

— Divide the space of the loop momenta into various re- 
gions and, in every region, expand the integrand in a 
Taylor series with respect to the parameters that are 
considered small there. 

— Integrate the integrand, expanded in the appropriate 
way in every region, over the whole integration domain 
of the loop momenta. 

— Set to zero any scaleless integral. 

To apply this strategy to a given limit one should first un- 
derstand, using various examples, which regions are rele- 
vant to it. In the Sudakov limit under consideration, these 
are the following regions, for a loop momentum k: 

hard (h): k ~ Q , 

M 2 

1-collinear (lc): fc+ ~ — — , fc_ ~ Q , k± ~ M , 



The basic identity of the Mellin-Barnes representation 
is the following (valid for | argAT — argF| < ir): 



2-collinear (2c): k + ~ Q , fc_ 

soft (s): fc~M, 
M 2 

ultrasoft (us): k < 



M . 



Q ■ 

By k ~ Q etc. we mean that any component of the vector k 
is of order Q, and k±, k± are the components of k denned 
after equation ijHJ. In other versions of the Sudakov limit, 
ultracollinear regions can also participate |84| . but they 
are irrelevant to the present version. 

So we obtain with this strategy the asymptotic expan- 
sion of our integrals as a sum of contributions generated 
by various regions. For brevity, we omit the word "gener- 
ated" and speak about contributions of regions, although 
the integration in each contribution is performed over the 
whole space of the loop momenta. 

In fact, this strategy is a generalization of the orig- 
inal strategy based on similar regions 44,39,41 , where 
the cut-offs specifying the regions were not removed so 
that the integrations were bounded by the regions under 
consideration. 



C Mellin-Barnes representation 

The Mellin-Barnes representation is a powerful tool for 
solving two closely related problems: a) The calculation of 
Feynman integrals, b) The asymptotic expansion of Feyn- 
man integrals in various kinematical limits. 



1 



1 



(X + Y) 



r(A) 



£n-^(A + ,)^. (108) 



It replaces a sum raised to any power by individual fac- 
tors which are raised to powers depending on the Mellin- 
Barnes parameter z. This simplification of the structure is 
obtained at the cost of an additional integration. The inte- 
gration contour in the Mellin-Barnes integrals runs from 
— ioo to -Moo and is chosen in such a way that poles from 
gamma functions of the form r(. . . + z) lie on the left hand 
side of the contour ("left poles") and poles from gamma 
functions of the form r(. . . — z) lie on the right hand side 
of the contour ( "right poles" ) . The contour cannot always 
be chosen as a straight line, especially if Re A < 0. 

If \X\ < \Y\, the integration contour can be closed on 
the left hand side at Kez — — oo, and the integral is given 
by the sum over the residues at the left poles. This sum 
corresponds to the Taylor expansion of 1/(X + Y) x for 
|X| < \Y\. On the other hand, if \X\ > \Y\, the contour 
can be closed on the right hand side at Rez = +oo, and 
the sum over the residues at the right poles corresponds 
to the Taylor expansion for |X| > \Y\. In the limiting case 
|X| = \Y\ the Mellin-Barnes integral is also convergent 
and is given by the sum of the residues on either the left 
or the right side of the contour - provided that these sums 
converge, which is often the case, especially when more 
than two gamma functions are present. 

The first application of the Mellin-Barnes represen- 
tation was, probably, in |SS]. The simplest possibility of 
using it is to transform massive propagators (X = k 2 , 
Y = —A/ 2 ) into massless ones (see e.g. ]86,87,88 as early 
references). In general, one starts from Feynman, alpha 
or Schwinger parameters and uses the Mellin-Barnes rep- 
resentation to separate arbitrary terms raised to some 
powers in such a way that the resulting parametric in- 
tegrals can be calculated in terms of gamma functions 
(see e.g. jSniEOlElES)- I n the context of dimensional reg- 
ularization, when the explicit evaluation at general values 
of d = 4 — 2s is hardly possible and one is oriented at 
calculating Feynman integrals in a Laurent expansion in 
e, the systematic evaluation by Mellin-Barnes represen- 
tations was initiated in |93II94| . An essential step of the 
evaluation procedure is the resolution of singularities in e, 
with the goal to represent a given multiple Mellin Barnes 
integral as a sum of integrals where the Laurent expan- 
sion of the integrands becomes possible. This is achieved 
by taking residues and shifting contours. Two different 
strategies for implementing this step were suggested in 
[9*5| and respectively. 

The identity l|108[) is valid for all powers A. In fact, 
the crucial point is not the convergence of the integral in 
the basic identity (|108|l . but the interchange of the order 
of integrations between the Mellin-Barnes integral and 
the (Feynman, alpha or Schwinger) parameter integrals. 
The necessary convergence of the parameter integrals re- 
stricts the real part of the Mellin-Barnes parameter z 
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to a specific range. If this range has a non-empty over- 
lap with the interval (— Re A, 0), the integration contour 
over z can be chosen as a straight line parallel to the imag- 
inary axis within the allowed range on the real axis. One 
can find values for the power A and other parameters such 
that an allowed range for the real part of z exists. The 
analytic continuation to the desired parameter values is 
then obtained by accounting for the residues which cross 
the fixed integration contour when the parameter values 
are smoothly changed .94,. Alternatively, the contour of 
the Mellin-Barnes integration can be deformed in such a 
way that it separates the poles of gamma functions with 
a dependence from the ones with a "— z" depen- 

dence even for the desired parameter values |Mj ■ Not only 
the gamma functions from the Mellin-Barnes representa- 
tion I|1U8|) , but also the ones introduced by the evaluation 
of the parameter integrals have to be considered here. As 
long as the prescription following equation (|108J) for the 
integration contour is respected, the convergence of the 
integrals is provided and all residues are accounted for on 
the correct side of the contour. This is still true if multiple 
Mellin-Barnes integrals are introduced by the iterated ap- 
plication of p()8|l . Even when A is a non-positive integer 
and -T(A) in the denominator gets singular, the right hand 
side of (|108|l is given by the limit where A approaches its 
actual value. In this case only a finite number of residues 
give non- vanishing contributions and reproduce the Bino- 
mial formula for (X + Y)I A L 

Often the Mellin-Barnes representation is used for 
asymptotic expansions (see e.g. |%Hll90H91 |l9"2 95.96 V 
When the Mellin-Barnes integral contains the factor t z 
with some parameter t, the asymptotic expansion in the 
limit t — > is given by picking up the residues on the right 
hand side of the integration contour. The asymptotic ex- 
pansion in the limit t — ► oo is given by the residues on 
the left hand side of the contour. By expanding t z in z 
about the poles of the integrand, the explicit form of the 
asymptotic expansion in powers of t and lnt with coeffi- 
cients from the Laurent expansion of the Mellin-Barnes 
integrand can easily be obtained In practice, how- 
ever, the most adequate way to perform the asymptotic 
expansion depends on the specific problem. For the work 
presented in this paper we have applied the method of ex- 
pansion by regions (appendix El and used Mellin-Barnes 
representations in the purpose of asymptotic expansion as 
a cross-check. When calculating scalar integrals for general 
propagator powers n, as in sectional the leading contribu- 
tions can be obtained from the Mellin-Barnes representa- 
tion by taking the residue at the first pole of each gamma 
function on the correct side of the integration contour. 
It turned out that in many cases the expressions obtained 
by the expansion of the loop integral within the expansion 
by regions method were simpler than the expressions ex- 
tracted from the Mellin-Barnes representation of the full 
integral. 

If the Mellin-Barnes representation is applied to the 
calculation of Feynman integrals (in particular, of indi- 
vidual contributions in an asymptotic expansion, as in 
the present work), when no large or small parameter t 



is present as t z in the Mellin-Barnes integrals or when 
the full dependence on t is desired, all residues on one 
side of the integration contour have to be considered and 
summed up. Some integrations in multiple Mellin-Barnes 
integrals can be performed explicitly by the application of 
identities based on the first Barnes lemma 97 , 

zoo 

/ ^ r ( A i + ^) r (A 2 + ^r(A 3 -z)r(A 4 -z) 

— ioo 

_ r(Ai + A 3 )r(Ai + A 4 )r(A 2 + A 3 )r(A 2 + a 4 ) 

r(A x + A 2 + A 3 + A 4 ) ' 1 j 

or on the second Barnes lemma |98| . 

zoo 

r dz r{\ 1 + z)r(\ 2 + z)r(\ 3 + z) 

J 2ni T(Ai + A 2 + A 3 + A 4 + A 5 + z) 

—too 

x r(A 4 - z)r(A 5 - z) 

_ r(Ai + A 4 )r(A 4 + A 5 )f(A 2 + A 4 )r(A 2 + A 5 ) 
r(Ai + A 2 + A 4 + A 5 )r(Ai + A 3 + A 4 + A 5 ) 
T(A 3 + A 4 )r(A 3 + As) 
r(A 2 + A 3 + A 4 + A 5 ) [1W) 

(see a collection of such formulae in Appendix D of |56|). 

Mellin-Barnes integrals develop singularities when a 
left pole and a right pole glue together in one point for 
some limit, e.g. e — > from dimensional regularization. 
These singularitities are directly present in the formulae 
(|109|) and (|110|l of the first and second Barnes lemma. In 
more complicated cases, it is usually a good idea to first 
extract the potentially singular residues by shifting the 
integration contours |93| or by an analytic continuation 
as described above and in Then the integrand may 
be expanded in the desired limits of its parameters. For an 
analytical result the residues on one side of the integration 
contour are summed up with the help of computer algebra 
programs, summation tables (see e.g. |56p or algorithms 
like [Bl)llrJrJ| . 

Characteristic examples of recent sophisticated calcu- 
lations based on the technique of Mellin-Barnes repre- 
sentations can be found in |99U1UU| . These results were 
crucial to check (in |1(J(J| ') cross order relations in N — 4 
supersymmetric Yang-Mills theory conjectured in |101j . 
Very recent results on checking the iteration structure in 
this theory with the help of Mellin-Barnes representations 
have been obtained in 102 103 104. 

Also recently algorithms for the automatic evaluation 
of Mellin-Barnes integrals have been formulated |lU5lTl 06 . 
These rely on the strategy of [IM] for the analytic continua- 
tion in the parameter e. The algorithms provide a basis for 
the analytic evaluation, and at least they can be applied, 
in their present form, to the numerical evaluation. The al- 
gorithm of jlUfcij is already implemented in Mathematica 
and would have been applied by us at least for numerical 
checks if it had existed early enough. 
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D Contributions in a theory with a mass gap 

The separation of the infrared-divergent electromagnetic 
contributions as described in 27 requires the two-loop 
corrections in the combined SU{2) x U(l) or U(l) x U(l) 
theory with massive and massless gauge bosons. In addi- 
tion to the results presented in the sections [3] and 0] of this 
paper, two-loop vertex and self-energy corrections with 
one massive SU{2) or U{\) gauge boson and one massless 
U(l) gauge boson are needed. 

As we regard an SU(2) x U(l) model without mixing 
between the two gauge groups (see pISlIM] for a discussion 
of this aspect), only the Abelian vertex and self-energy 
diagrams ffiguresEland ll lfl contribute. After replacing one 
of the two massive SU{2) gauge bosons in these diagrams 
by a massless U(l) gauge boson, we obtain the results 
listed in the following paragraphs. 

The planar vertex correction of the diagram in fig- 
ure^) with line 5 (cf. figure massless is 



J - 1 M 5 =0 _ n 




4£ 2 



3 



2 2 A , 7 2 85 



+ -£ 4 - -£ 3 
3 3 



17 )£ 2 



2 2 101\ „ 107 4 184 

-36C 3 + -7T 2 ) C H 7T 4 H C 3 

S3 3 2 / 90 3 s3 



59 2 599 m . . ^(M 2 
T2* + ~ \ +0{£)+ °W 



(111) 



When either of the two gauge bosons in the non-planar 
vertex diagram of figure[2b) is massless (cf. figurc[S]for the 
line numbering), the contribution is 



P M 5 =0 
f v,NP 



P M 6 =0 
r v,NP 



2 \ 2e 

V-\ s 2 

(4tt) 2 U 2 1 ' 



--£ 3 + AC 2 - 12C - 12Cs + 7T 2 



lc 4 - 4£ 3 



3 



22 \C 2 



^ 2 -68)c-^ 
3 / 90 



56Cs 

4tt 2 + 96 \ + 0(e) + O 



(Ml 
\Q 2 



14 



90Cs 
(113) 



The Mercedes-Benz graph in figure |2t) gives the contri- 
bution 



",BE - C F 




89\ 13 4 

— £ 7r + 16Ca 

2 90 ^ 3 



, (H4) 



when line 3 (cf. figure EJ) is massless, and 



where a and a' are the couplings of the SU{2) and U(l) 
gauge groups, respectively, and £ = ln(Q 2 /M 2 ). The 1/e 3 
pole is due to the infrared divergence. The same diagram 
with line 6 massless yields the contribution 



i-,M 6 =0 _ n 



(4tt) 2 \M 2 
1 




3 
13 



-£ 2 + 3£ - 



1 £ 2 



11 

~ T 

-24C 3 - ^ 2 



:C 4 



11 

T 



4 13 2 41 

H 7T 4 + 46C 3 + — 7T 2 

45 S 12 8 



/ M 2 



(112) 



Note that only the linear logarithm and the non-logarith- 
mic constant at order e° of this result differ from the 
case (11411 with two massive gauge bosons (and, of course, 
the different prefactor Cf act' instead of C F a 2 ). 



jpM 4 =0 _ n 



EL 

(4tt) 2 \M 2 
1 





2 1 










e 3 e 2 


2. 



-2 , 7T 53 

-It 1 + 7C- — - — 



11 4 32 ^ 17 2 479 

H 7T 4 C3 + — 

45 3 S 4 8 



f2 



+ 0(s) + 0, 



(115) 



when line 4 is massless. 

The vertex correction with fermion self-energy of the 
diagram in figure |5Ji) yields 



F, 



M 3 =0 



v,fc 



= C F 



aa 

(4^)2 VM 2 J 
1 



2= 



1 



-2£ 



2 £ 2 _ 7 £+^ + ^ 



- -£ 3 + 7£ 2 
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\Q 2 



(116) 



with line 3 (cf. figure 0) massless and 

2e ( 



r M 5 =0 
v.fc 



aa! f ^ 2 V 



(4tt) 2 

- - L- ■ - + 
e [ 3 

r 2 



1 



£ 2 - 3£ + V + - 

4 



-£ 3 



5£ 2 



7T 

4C3+ 12 



163 



0{e)+0 



M 2 

Q 2 



(117) 



with line 5 massless. The only difference of (|117fl with 
respect to the purely massive result (|49|) is in the non- 
logarithmic constant at order e . 

The self-energy diagram in figure ITTk ) contributes 



■"vx 



Tl 



C f 



(4tt) 2 
-0{e) 



i^V% 2 (L-l 

M 2 £ \2e 4 



(118) 



when either of its two gauge bosons is massless. The two 
contributions of the self-energy diagram in figure ITTb ) are 



-tM2=Q 

J T2 



= C f 



(4tt) 2 
-0(e) 



M 2 I E \ 2e 2 4e 4 8 



(119) 



with the gauge boson in the outer loop massless, and 



Zj T2 



C, 



aa 
(4^)2 

■0(e) 



M 2 



2e 



2e 2 



1 
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T 8, 
(120) 



with the gauge boson in the inner loop massless. Note that 
112UI) differs from the purely massive result l|76|) only at 
order e° . 

The self-energy diagram in figure II lb ) has no corre- 
sponding contribution with one massive and one massless 
gauge boson, because these integrals vanish in dimensional 
regularization. 

According to (|69(l , the product of the massless one- loop 
vertex correction F^\ =0 and the (massive) one-loop self- 
energy correction Si l|71|l is needed as well. The missing 
piece is well known: 




+ 0{e 2 ). (121) 



The prefactor has to be expanded as 



7TT) =(^) [)-:C + -f--C>) + 0(:> } 



in order to match the prefactor of the other contributions. 

The contributions to the SU(2) x U(l) form factor 
with one massive and one massless gauge boson may now 
be added together in analogy with Q99|l : 



Fi 



?M 5 =0 



pM 6 =0 
r v,LA 



M 5 =0 



P M 6 =0 
r v,NP 



= c 




M 5 =0 
v,fc 



12C 2 



(122) 



The infrared-convergent two-loop interference term of 
equation (6) in results from l|122|l after subtraction 
of the massive times the massless one-loop form factor: 
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